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Abstract _ _ _ 

This  report  defines  a  dimensionless  shape  factor,  which  is  useful  for  characterizing  the 
penetration  potential  of  projectiles  and  irregular  shaped  fragments.  The  shape  factor,  so 
defined,  is  purely  a  function  of  shape  and  orientation  and  is  independent  of  mass  and  material 
density.  The  shape  factor  of  some  simple  shapes  is  calculated  exactly,  but  the  focus  of  this 
report  is  on  a  right  circular  cylinder  (RCC).  Not  only  is  it  possible  to  express  the  shape 
factor  as  a  function  of  the  length-to-diameter  (L/D)  ratio  and  orientation,  but  it  is  also 
possible  to  derive  an  exact,  closed-form  expression  for  the  shape  factor  probability 
distribution.  It  is  found  that  the  probability  density  function  is  not  a  symmetrical  distribution 
about  its  mode,  but  rather  is  highly  skewed.  This  points  out  the  inadequacy  of  an  average 
shape  factor  and  also  carries  implications  for  designing  fragment  simulating  projectiles 
(FSPs).  Furthermore,  it  is  shown  that  randomly  oriented  cylinders  have  potential  for 
simulating  behind-armor  debris  fragments. 
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1.  INTRODUCTION 

This  report  defines  a  dimensionless  shape  factor,  which  is  useful  for  characterizing  the  penetration  potential 
of  projectiles  and  irregular  shaped  fragments.  The  shape  factor,  so  defined,  is  purely  a  function  of  shape  and  orienta¬ 
tion  and  is  independent  of  mass  and  material  density.  This  is  a  departure  from  the  past  where  the  previous  definition 
incorporated  the  material  density— so  that  the  shape  factor  for  a  steel  fragment  was  different  from  an  aluminum 
fragment,  even  when  they  had  the  same  size  and  shape.  Indeed,  this  fact  is  easy  to  overlook  with  the  past  definition, 
and  so  one  can  get  erroneous  results  by  applying  the  shape  factor  for  one  fragment  material  to  another  fragment 
made  from  a  different  material.  Past  treatments  of  shape  factor  have  been  inadequate  in  two  other  respects: 

(1)  Average  shape  factor  (i.e.,  the  mean  value  over  all  orientations)  has  been  used  to  characterize  the  entire 
shape  factor  distribution  over  random  orientations. 

(2)  To  improve  upon  the  average  shape  factor,  it  is  common  practice  to  use  a  normal  distribution  to  simulate  an 
arbitrary  orientation. 

This  report  challenges  these  two  assumptions.  We  show  that  it  is  possible  to  calculate  the  exact,  closed-form  expres¬ 
sion  for  the  shape  factor  probability  distribution  over  all  orientations  of  a  right  circular  cylinder  with  arbitrary  length 
to  diameter  ratio.  This  reveals  that  average  shape  factor  does  not  do  a  good  job  of  characterizing  an  arbitrary  orien¬ 
tation,  and  that  the  probability  distribution  is  neither  normal  nor  symmetric  but  rather  is  skewed  in  favor  of  larger 
shape  factors.  Since  fragment  shape  factor  plays  such  an  important  role  in  penetration  and  perforation,  an  improved 
treatment  of  shape  factor  will  have  important  consequences  for  predictions  of  fragment  perforation  and  for  the 
design  of  fragment  simulating  projectiles  (FSPs). 

Fragment  presented  area  upon  impact  with  a  target  plate  is  a  key  determinant  for  penetration,  as  is  fragment 
mass.  Indeed,  apart  from  velocity,  the  most  important  factor  for  penetration  is  probably  the  mass  per  unit  presented 
area.  This  is  borne  out  by  an  examination  of  the  THOR  penetration  equations  [1,2].  The  THOR  velocity  retarda¬ 
tion  equation  is 

vr  =  vs  -  Icv(tA)am f  (sec  6)yvxs, 

where  vr  is  the  residual  velocity,  v*  is  the  striking  velocity,  kv  is  a  constant,  t  is  the  target  plate  thickness,  A  is  the 
average  presented  area  of  the  projectile,  ms  is  the  striking  mass,  and  9  is  the  obliquity  angle.  The  exponents  a  and 
P  depend  upon  the  target  material.  Figure  1  is  a  plot  of  these  exponents  for  20  different  target  materials,  ranging 
from  very  soft  to  very  hard. 
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Figure  1.  Correlation  Between  Mass  and  Presented  Area  Exponents  in  THOR  Penetration. 

It  is  evident  that  experimental  data  support  the  notion  that  —p  =  a,  according  to  the  dashed  line — which 
means  that  mass  and  presented  area  only  enter  into  the  THOR  penetration  equation  in  the  specific  combination  of 
mJA.  Projectile  mass  is  universally  recognized  as  a  key  parameter  for  penetration,  whereas  presented  area  is  per¬ 
haps  not  appreciated  to  the  same  extent.  But  we  see  here  that  they  are  equally  important  in  the  THOR  penetration 
equation. 
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Aside  from  its  volume  (which  is  determined  by  mass  and  density),  the  presented  area  of  a  projectile  is  com¬ 
pletely  characterized  in  terms  of  shape  and  (except  for  spheres)  orientation.  Both  of  these  are  contained  in  the 
notion  of  shape  factor,  which  we  will  define  momentarily.  In  this  report,  we  restrict  ourselves  to  the  shape  of  right 
circular  cylinders,  for  three  primary  reasons: 

(1)  Fragment  simulating  projectiles  are  usually  cylindrical  for  ease  of  firing  from  a  gun  barrel  (as  well  as  the 
ability  to  control  the  impact  area  through  spin-stablization). 

(2)  Fragment  penetration  models,  such  as  FATEPEN  (Fast  Air  Target  Encounter  Penetration)  [3],  often  use 
cylinders  to  simulate  the  fragment  shape. 

(3)  It  is  simple  enough  to  possess  an  exact,  closed-form  expression  for  the  probability  density,  yet  complex 
enough  to  have  potential  for  simulating  behind-armor  debris  fragments. 

This  report  focuses  on  right  circular  cylinders  with  arbitrary  Li D  ratios,  where  L  is  the  cylinder  length  and  D  is  the 
diameter.  Also,  both  specific  orientation  (e.g.,  yaw  angle)  as  well  as  random  orientation  (uniform  over  the  surface  of 
the  unit  sphere)  are  treated. 

2.  DIMENSIONLESS  SHAPE  FACTOR 

The  presented  area,  Ap,  is  a  function  of  the  dimensions  of  the  cylinder,  its  shape,  and  its  orientation.  The 
only  combination  of  the  mass  and  material  density  that  has  the  dimensions  of  an  area  is  (mass/density)  .  There¬ 
fore,  the  presented  area  is  of  the  form 

Ap  =  y(m/p)2ri,  (!) 

where  y  is  a  dimensionless  shape  factor  that  accounts  for  both  fragment  shape  and  orientation.  It  is  purely  a  func¬ 
tion  of  geometry  and  is  independent  of  both  mass  and  density. 

Another  definition  of  fragment  shape  factor  that  has  been  used  rather  extensively  is  defined  by  the  equation 

Ap  =  s  /n273.  .  (2) 

To  contrast  this  with  the  dimensionless  shape  factor,  we  call  this  quantity  s  the  conventional  shape  factor.  However, 
this  is  not  a  good  way  to  define  a  shape  factor.  Aside  from  the  fact  that  the  units  are  the  rather  cumbersome 
cm2/gm2'3,  or  in^grains273,  the  shape  factor  so  defined  is  dependent  upon  the  material  density.  Thus,  if  we  deter¬ 
mine  the  shape  factor  for  steel  fragments,  then  they  cannot  be  applied  to  tungsten  fragments  even  if  they  have  the 
same  shape  without  first  accounting  for  the  difference  in  densities.  The  dimensionless  shape  factor  defined  by  equa¬ 
tion  (1)  does  not  have  this  restriction.  It  is  easy  to  convert  old  values  to  dimensionless  values  since  the  relation 
between  the  two  is  simply 

[  p273  s  if  s  is  in  units  of  cm2/gram2/3 
^  1  40 p273  s  if  s  is  in  units  of  in2/grains273 

where  1  gram  =  15. 432  grains.  Since  the  dimensionless  shape  factor  is  purely  a  function  of  geometry,  it  is  possible 
to  calculate  it  for  various  standard  shapes,  as  shown  in  Table  1 . 
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Table  1.  Dimensionless  Shape  Factors  for  Some  Common  Shapes 


Shape 

Orientation 

Shape  Factor ;  y 

Sphere 

Not  applicable 

(3/2)2l3(!T/4)m  ~  1.209 

Cube 

Face-forward 

1 

Cube 

Edge-forward 

t/2®  1.414 

Cube 

Random3 

3/2 

Cube 

Comer-forward 

V3-  1.732 

Cylinder  (L/D=l) 

Face-forward 

(^)1/3~  0.927 

Cylinder  (L/D=l) 

Side-forward 

(5)-2/3  *  1.175 

Cylinder  (L/D=l) 

Random3 

(3/2)(tt/4)  1/3 «  1.384 

Artillery  Fragment15 

Random3 

2.1 

Spall  Fragment6 

Random3 

2.3 

a  For  random  orientation,  the  rule  Ap  =  (Surface  Area)/4  was  used  (see  Appendix  A), 
b  The  values  for  artillery  and  spall  fragments  were  derived  from  the  measured  presented  area. 


In  the  remainder  of  this  report,  we  use  the  terminology  “shape  factor”  to  mean  the  dimensionless  shape  factor,  as 
defined  by  equation  (1). 

3.  SHAPE  FACTOR  OF  A  CYLINDER 

The  orientation  of  the  cylinder  can  be  parametrized  by  the  angle  0 ,  as  shown  in  Figure  2. 

I  Axis  of  Symmetry 


Viewing  Direction 


Figure  2.  Right  Circular  Cylinder  With  Orientation  Angle, 

When  viewed  from  the  side  (9  =  ni 2),  the  projected  area  is  LD.  When  viewed  from  the  top  (9  =  0),  the  projected 
area  is  kD21  4.  For  an  arbitrary  direction,  9 ,  the  projected  area  is  the  linear  combination  of  these  two  faces  scaled  by 
the  cosine  of  the  projected  area  surface  normal.  Thus, 

Ap  =  LD  cos(x/2  -  0)  +  j  D2 1  cos  0!  =  D2  [L  sin  0  +  ^  I  cos  01  j.  (4) 

It  is  necessary  to  take  the  absolute  value  of  the  cosine  term  because  cos  0  is  negative  when  n/2  <9  <k.  The  volume 
of  the  cylinder  is  given  by 

ml  p  =  1  D2L.  (5) 

Substituting  these  two  expressions  into  equation  (1)  and  solving  for  y,  we  find  that 
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(6) 


(n  L\  (L  it  A 

r  \AD)  U>  4  J 

It  is  convenient  to  define  the  parameters 


Then  the  shape  factor  may  be  written  as 

y(0)  =  a  sin  9  +  b\  cos  0\, 


(7) 

(8) 


where  0  <9  <n.  Notice  that  the  presented  area  of  the  cylinder  has  the  same  value  at  the  viewing  angle  9  as  it  has  at 
K  _  fft  so  that  y(9)  =  yin  -  9).  Thus,  we  could  restrict  the  orientation  to  0  <  9  <  nil  and  just  perform  averages  over 
the  top  hemisphere  of  the  unit  sphere.  This  would  allow  us  to  drop  the  absolute  value  signs  in  equation  (8).  Never¬ 
theless,  to  avoid  confusion,  we  shall  use  the  form  expressed  in  equation  (8)  and  perform  averages  over  the  entire  unit 
sphere. 

This  expresses  the  shape  factor  as  an  explicit  function  of  the  orientation  angle,  9.  It  is  also  an  implicit  func¬ 
tion  of  the  LID  ratio  through  the  coefficients  a  and  b.  A  plot  of  7  as  a  function  of  both  orientation  angle  9  and  LID 
ratio  is  shown  in  Figure  3. 


Figure  3.  Shape  Factor  of  a  Cylinder  as  a  Function  of  Orientation  Angle  and  LID  Ratio. 


From  equation  (6),  we  see  that  7(0)  =  b  and  y(nll)  =  a.  The  minimum  value  of  7  is  one  of  these,  depending  upon 
the  particular  LID  ratio: 

f  a  if  LID  <  nl4 


7min  —  ‘ 


b  if  LID  >  n/4 


or,  simply,  ymin  =  min(a,  b). 


The  maximum  value  of  7  can  be  obtained  by  setting  dy!d9  —  0,  and  this  gives 


7max  =  Va2  +  &2  =  (^-| 


-2/3! 


K 

-|+ID 


-|l/2 


The  maximum  shape  factor  is  realized  at  the  orientation  angle 


*r=jw=tan 


if  W4"£l 


(9) 


(10) 


(11) 


Thus,  in  general,  there  are  two  graphs  for  the  shape  factor  as  a  function  of  orientation  angle,  depending  upon  the 
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LID  ratio  of  the  cylinder.  These  two  cases  are  exemplified  in  Figures  4  and  5. 


Figure  4.  Shape  Factor  Dependence  Upon  Orientation  for  LID  >  id  A. 


Figure  5.  Shape  Factor  Dependence  Upon  Orientation  for  LID  <  7tfA. 

The  angle  0y=Q  in  Figure  4  and  0Y=b  in  Figure  5  can  be  obtained  through  the  use  of  equation  (8)  and  we  find 

07=a  =  tari ~l{a  J?  )  when  LID  >  id  A  (12) 

and 

0y=b  =  tan-1(  -- )  when  LIDkttIA.  (13) 

—  a1 

From  Figures  4  and  5,  we  also  see  that  when  ma x(a,  b)  <  y  <  7max,  the  orientation  angle  is  double-valued,  and  when 
min  (a,  b)  <  7  <  max(a,  b),  the  orientation  angle  is  single-valued. 

4.  SHAPE  FACTOR  AS  A  FUNCTION  OF  ORIENTATION 

The  minimum  shape  factor  is  given  by  equation  (9),  and  the  maximum  shape  factor  is  given  by  equation 
(10).  The  mean  shape  factor,  when  averaged  over  all  orientations  equally,  is  given  by 
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In  n 


=  J  y{8)dQ.  =  J  d<t>  J  y{0)  sin  0 dd. 


(14) 


Substituting  equation  (8)  and  performing  the  integration,  we  find+ 

n  1  f/rY'Yl 


■2/3 


L  1 
D  +  2 


(15) 


Calculation  shows  that  y  is  a  minimum  when  LID  =  1,  where  it  takes  on  the  value  (3/2) (^/4) 1/3  =  1.38395.  If  we 
only  know  the  average  shape  factor  of  a  fragment,  and  it  exceeds  this  value,  then  eq.  (15)  can  be  used  to  find  the  LID 
ratio  of  an  equivalent  disk-shaped  cylinder.  This  solution  is  shown  plotted  in  Figure  6. 


Figure  6.  Disk-Shaped  Cylinder  Derived  From  the  Average  Shape  Factor. 


Here  is  a  simple  program  that  implements  this  solution  of  the  cubic  equation  for  the  LID  ratio: 


//  asf21d.C:  Find  the  L/D  ratio  of  an  ECC  with  the  given  average  shape  factor. 

#include  <iostream.h> 

#include  <math.h> 

void  main<  void  ) 

const  double  C  —  pow(  0.25  *  M_PIf  -1.  /  3.  ); 

const  double  G_BAR_MIN  =  1.5  /  C; 

double  l_dMin,  gBar; 

while  (  cin  >>  gBar  )  { 

if  (  gBar  >  G_BAR_MIN  )  [ 

double  p  =  sqrt(  2 .  *  C  *  gBar  /  3 .  )  ; 

double  q  =  acos(  ~pow(  pf  -3.  )  ); 

l__dMin  =  pow (  2.  *  p  *  cos(  q  /  3.  ),  -3.  )? 

else  if  (  gBar  ==  G_BAR_MIN  )  l_dMin  =  1. ; 
else  { 

cerr  <<  "The  minimum  average  shape  factor  for  an  RCC  is 
<<  G_BAR_MIN  <<  endl ; 
exit(  1  ); 

} 

cout  <<  l_dMin  <<  endl; 

) 

} 


This  result  also  follows  from  the  general  theorem  that  the  average  projected  area  of  a  convex  solid  is  one-fourth  the  total  surface  area  (see 
Appendix  A). 
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The  average,  minimum,  and  maximum  shape  factors  are  plotted  in  Figure  7. 


0  2  4  6  8  10 

LID 

Figure  7.  Minimum,  Maximum,  and  Mean  Shape  Factors  as  a  Function  of  LID  Ratio. 

It  is  straightforward  to  establish  the  following: 

•  y  is  a  minimum  when  LID  -  1 . 

•  ymax  is  a  minimum  when  LID  =  V2  id  A  «  1.111. 

•  7min  is  a  maximum  when  LID  =  ar/4  «  0. 785. 

•  ymax  -  ym[n  (i.e.,  the  range  of  y  over  all  possible  orientations)  is  a  minimum  when  LID  =  ar/4  »  0. 785. 

Its  value  at  the  minimum  is  (V2-  l)(^r/4)^1/3  *  0.449.  This  quantity  is  a  measure  of  the  deviation  from 
spherical  symmetry,  since  a  sphere  has  ymax  -  ymin  =  0. 

For  completeness,  we  can  also  calculate  the  variance: 


nil 

a2  =  —  |  (/(6>)  -  ffdQ.  =  J  (y(6»)  -  f  )2  sin  0  dO, 


and  we  find 


2  2  f2  » 1  ,  1  ,2 

a  = - -  a  + - -  ab  +  —  b  . 

,3  16  J  1 3  4  J  2 


5.  SHAPE  FACTOR  DISTRIBUTION  OF  A  RANDOMLY  ORIENTED  CYLINDER 

It  is  possible  to  derive  a  closed-form  expression  for  the  shape  factor  probability  distribution  function.  We 
start  with  the  formula  for  the  probability  density. 


f(y)  =  J  S(asin0  +  b cosO -y)sin6d0y 


where  S(x)  is  the  Dirac  delta  function,  or,  setting  x  =  cos  0 , 


f(7)=  j  8(g{x))dx, 


where 


g(x)  s  aV  1  -  x2  +  bx  -  7.  (20) 

To  evaluate  this  integral,  we  use  the  fact  that  the  delta  function  of  a  function  can  be  expressed  in  terms  of  the  roots  of 
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the  function  [4]: 


where  xt  (i  =  !,•*■ 


n 


<*(£(*»  =  X 

1=1 

,  n)  are  the  n  roots  of  the  equation  g(x) 


S(X-  X;) 

I dgldx\x=x.  ’ 

=  0.  The  roots  of  g(x)  are  found  to  be 


by  ±  +  — 

a2  +  b2 


This  can  be  simplified  by  introducing  the  quantities 


7mm 


- ,  and  y‘ 

7mm 


7 

7m  ax 


where  ymm  is  given  by  eq.  (10).  The  roots  can  now  be  written  as 

x  -  b'y'  ±  a'Vl  -y'2. 

Now  notice  that  a’2  +  b'2  =  1  and  0  <  y'  <  1  so  that  we  can  define  angles 

a  =  cos-1(6')  =  cos-* (b/ ymax )  and  P  =  cos-I(y')  =  cos -1(y/ymax). 
This  allows  us  to  express  eq.  (24)  as 


x  =  cos  a  cos  P  ±  sin  a  sin  p  =  cos(a  ±  P). 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 


So  the  two  roots  are 

xi  =  cos(«  -  p)  and  x2  =  cos  (a  +  p), 

or,  in  terms  of  the  orientation  angle,  0, 

6>i  =  a-  p  and  02  =  a  +  p. 

Evaluation  of  the  derivative  of  g(x)  is  straightforward,  and  we  find 

,,  s  7mm  Sin(-P)  rmaxSin^ 

g  (x i)  =  — TT77 — —  311(1  8  (*2)  = 


sin(a  -  P) 


Thus, 


sin(a  +  P) 

-x{)  8{x~. 


(27) 

(28) 

(29) 


(30) 


and,  using  equation  (29), 


f(7)  = 


I  sin(a  -  P)  I 

l7max  Sin(-/?)l 


<5(x-X[)  + 


I  sin(a  +  p)  I 
l}W  sin(/5)l 


S(x  -  x2 ) 


(31) 


In  order  to  evaluate  this  integral,  it  is  sufficient  to  know  the  location  of  the  two  roots,  xx  and  x2,  or,  equivalently,  8\ 
and  02.  For  LID  >  n/A,  we  find  that  when  b<y  <  a,  then  0  <  0X  <  0y=a  and  when  a  <  y  <  ymax,  then  6>,  =  0y=a  and 
0r=Km  <02^  nil.  Similarly,  when  LID  <  n/A,  we  find  that  when  a<y  <b,  then  8y=b  <02<  nil  and  when 
b  <  y°<  ymax,  then  0  <  9X  <  0r=yma  and  02  =  0y=b.  Thus,  only  one  root  lies  in  the  interval  where  the  orientation 
angle  is  single-valued  and  both  roots  lie  in  the  interval  where  the  orientation  angle  is  double-valued.  The  evaluation 
of  the  integral  is  now  trivial,  and  the  results  are  as  follows. 


•  L/D>n/4 


If  b  <  y  <  a,  then  f(y)  = 


sin(a  -  P) 
7mm  sin  p  ' 


If  a  <  y  <  7max>  then  f(y)  = 


sin(a  -  p)  sin(a  +  P) 
7mm  sin  P  7mm  P 


(32) 

(33) 
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.  L/D<n/4 


If  a<y  <  b,  then  /(/)  - 


sin  (a  +  p) 
/max  sin  P 


If  b  <  7  <  7max>  Aen  /(/)  = 


sin(a  -  p)  sin(a  +  p) 

/max  Sin  P  /max  Sin  P 


Expressing  these  results  in  terms  of  the  parameters  a  and  b  gives  us  the  following  results. 


5.1  Probability  Density  Function 


.  L I D  >  n!  4 


/(r)  = 


«/-Wr^  ax-r2 

rmaxV^ax-72 

2a/ 


7maxV7max-72 
Figure  8  is  a  plot  of  this  function  for  a  compact  cylinder. 


b  <y  <  a  0 <6  <64 
a<y  <  /max  Bd<e<Jt!  2 


(34) 

(35) 


(36) 


0.9  1  1.1  1.2  1.3  1.4  1.5 


7 

Figure  8.  Shape  Factor  Probability  Density  Function  for  a  Cylinder  With  LID  >  k! 4. 

Notice  the  discontinuity  at  /  =  (tr/4)-273,  where  the  shape  factor  changes  from  being  single-valued  to  being  double¬ 
valued. 

.  L! D  <k! 4 


f(7)  = 


Q7  +  W7max~72 
7maxV7max  —  7 
lay 

7maxV7max-72 


a<y  <b  6d  <6  <  nil 
fc</</max  0<6<6d 


(37) 


Figure  9  shows  a  plot  of  these  functions  for  a  disk-shaped  cylinder. 
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Figure  9.  Shape  Factor  Probability  Density  Function  for  a  Cylinder  with  LID  <  jt/4. 

Notice  that  the  most  probable  value  of  the  shape  factor,  (i.e.,  its  mode),  is  ymax,  regardless  of  the  magnitude  of  the 
LID  ratio. 

5.2  Cumulative  Distribution  Function 


UD>it!4 


i _ frr  +  W/max  if  b<y< a  o<e<ed 

/max 

l _ 2W/max -_zL  a < 7 < ymax  ed<e<7tn 


Figure  10  is  a  plot  of  the  cumulative  distribution  for  a  compact  cylinder. 


i  H  F(r) 


0.9  1  1.1  1.2  1.3  1.4  1.5 


Figure  10.  Shape  Factor  Cumulative  Distribution  Function  for  a  Cylinder  With  LID  >  nt 4. 


•  Case:  L/D  <  7r  /  4 


F{y)  =  \ 


by-a-fy^-y2 


Jraax 


^V^max-y2 

/max 


a  <y  <  b  6(l<9<  k!2 


b<y<  ymax  o<e<ed 


A  plot  of  this  cumulative  distribution  for  a  disk-shaped  cylinder  is  shown  in  Figure  1 1 . 


(39) 


Figure  11.  Shape  Factor  Cumulative  Distribution  Function  for  a  Cylinder  With  LID  <  x/4. 


5.3  Generating  the  Random  Shape  Factor  Distribution 

The  cumulative  distribution  expressed  by  equations  (38)  and  (39)  can  be  inverted  to  give  y  in  terms  of  the 
cumulative  probability.  Then  by  selecting  the  cumulative  probability  from  a  uniform  random  distribution  on  the  unit 
interval,  we  generate  a  random  distribution  of  shape  factors,  where  the  frequency  of  occurrence  follows  the  probabil¬ 
ity  density  function — a  standard  technique  for  generating  random  numbers. 

•  Case:  UD >k/ 4 


b(\  —  P)  +  a'JT~—{l~—-py2  0<P<\-2ab/y- 


2 

max 


(40) 


ymaxV l-(ymax(l--P)/2fl)2  1  -  2ab/y^  <  P  <  1 


•  Case:  L/D  <  x  /  4 


y  = 


bP  +  crfl  -P2 


0  <P<(b2-  a2)/yt 


2\/,,2 

max 


ymaxVl-(ymaxO-^)/2a)2  (b2  -  d2)^  <  P  <  1 


(41) 


These  equations  enable  us  to  define  an  algorithm  for  generating  random  shape  factors  having  the  same  distribution 
as  that  corresponding  to  the  uniform  random  orientation  of  a  cylinder.  First,  define  the  following  quantities: 
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7t  L  JU  7t  Lt  Tt  j  «  "  7*9*  i 

-  —  b=  --  ymax  =  ^a2  +  b2,  P ,sl- 


4  Z> 


r2/3 


D 


k  L 


4  D 


•2/3 


2tf£ 

/max 


and  P2 


(fc2 


a2) 


/max 


Then  the  algorithm  is: 

(1)  Generate  P  ~U (0,1). 

(2)  If  a>b,  (i.e,  LID  >  7t!A), 

(2a)  If  P  <  Pi ,  return  y  =  b(l  -  P)  +  flVl  -  (1  ~?)2- 
(2b) Else  if  P>Pi9  return  y  =  /maxV  1  ~  (ZmaxO  -  P)l2a)\ 

(3)  Else  if  a  <  b,  (i.e,  LID  <  k!A\ 

(3a)  If  P  <  P2,  return  y  =  bP  +  aV  1  -  P2.  _ 

(3b) Else  if  P  >  P2,  return  y  =  /maxV1  -(/maxO  -P)^a)2. 

This  algorithm  is  more  than  twice  as  fast  in  execution  time  than  the  considerably  simpler  shape  factor  simulation : 

(1)  Generate  ( 9 ,  <p)  ~  Surface  of  the  Unit  Sphere.1 

(2)  Return y  =  asin6  +  b cos 0. 

Before  moving  on,  we  summarize  some  properties  of  the  shape  factor  distribution  for  a  randomly  oriented  cylinder 
in  Table  2. 


Table  2.  Shape  Factor  Parameters  for  a  Randomly  Oriented  Cylinder 


Parameter 


Expression 


Values  when  LID  =  1 


Length 

Diameter 

Parameter  a 

Parameter  b 

Minimum 

Maximum 

Mean 

Variance 

Mode 


Median 


L 

D 


t  v-2/3  j 
n  L  I  L 


*  =  14  D 


D 


i  \-2/3  _ 
k  L\  n 


b  =  '4D)  4 

7min  =  min  (a,  b) 

/max  =Va2  +  fc2 

-  K  1  u 
y  =  —  a  +  -b 

4  2 


^  l  2 

1^3  16  J  [3  4  J  12 


/  “  /max 


/  = 


/maxV1  ~(7mJ4a)2  if  4  <  p  <  ^ +  2)  4 

if(V5  +  2)f<§ 

/maxV^  —  ^  ^  ~4  <  5  <  ~4 

5ia  +  b  L  1  n 

-  if  —  <  — 

2  £>  V3  4 


L 

L 

a  =  l.  17474 

Z>  =  0.922635 

7min  =  0.922635 
7maj  =  1.49374 

7  =  1.38395 

a2  =  0.0109974 
7  =  1.49374 


7  =  1.41626 


t  See  Saucier  [5]  for  an  algorithm  to  generate  angle  pairs  (8, 6)  that  lie  on  the  surface  of  the  unit  sphere. 
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6.  SIMULATION  OF  BEHIND-ARMOR  DEBRIS  FRAGMENTS 

As  the  LID  ratio  decreases,  and  the  fragment  becomes  more  disk-like,  the  probability  density  function  f(y) 
becomes  nearly  constant  everywhere  except  very  close  to  ymax,  where  it  becomes  highly  peaked.  Consequently,  the 
cumulative  distribution  function  F(y)  becomes  almost  a  straight  line  over  nearly  its  entire  domain.  For  example, 
Figure  12  shows  a  plot  of  the  F(y)  for  an  LID  of  0. 1 . 


Figure  12.  Shape  Factor  Cumulative  Distribution  Function  for  a  Cylinder  With  Small  LID . 

Thus,  if  fragments  can  be  simulated  with  small-L/D-ratio,  disk-shaped  cylinders,  then  we  should  expect  a  relatively 
flat,  or  uniform,  shape  factor  distribution  over  orientation  angle.  While  this  pertains  to  a  single  fragment,  what  about 
a  distribution  of  fragments?  Since  there  is  no  obvious  reason  why  the  fragments  would  be  expected  to  have  all  the 
same  LID  ratio,  let  us  assume  that  they  vary  over  some  finite  range.  How  will  this  affect  the  shape  factor  distribu¬ 
tion?  For  example.  Figure  13  shows  four  cylinders  with  different  LID  ratios.  As  a  point  of  reference,  a  dime  has  an 
LID  of  0.1. 


LID  =  1 


LID  =0.5 


LID  =  0.1  LID  =  0.05 


Figure  13.  Four  Cylinders  of  the  Same  Mass  but  Different  LID  Ratios. 
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Let  us  assume  that  the  LID  ratio  is  equally  likely  to  be  anywhere  from  0.05  to  0.5,  and  that  each  fragment  is  oriented 
uniformly  over  the  unit  sphere.  To  determine  the  resultant  shape  factor  distribution,  we  perform  a  Monte  Carlo  sim¬ 
ulation.  Figure  14  shows  a  typical  result. 
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Figure  14.  Shape  Factor  Simulation  from  Randomly  Oriented  Cylinders  With  Uniformly  Distributed  LID  Ratios. 


The  distribution  appears  to  be  lognormal.  In  any  case,  it  is  skewed  to  larger  shape  factors,  and  is  far  from  uniform. 
In  summary: 

•  A  fragment  in  the  shape  of  a  disk-like  cylinder  with  a  small  but  fixed  LI D  ratio  leads  to  a  uniform  distribu¬ 
tion  of  shape  factors. 

•  A  collection  of  fragments  in  the  shape  of  disk-like  cylinders  with  small  but  varying  LID  ratios  leads  to  a 
distribution  of  shape  factors  that  is  skewed  to  larger  values  and  appears  to  be  lognormal. 

The  next  step  is  to  compare  these  results  to  experimental  data. 

6.1  Experimental  Design  to  Collect  Behind-Armor  Debris 

Behind-armor  debris  fragments  were  collected  in  Cellotex  when  a  30-mm  armor  piercing  discarding  sabot 
(APDS)  round  was  shot  at  a  1-inch-thick  target  of  rolled  homogeneous  armor  (RHA).  The  debris  fragments  were 
subsequently  collected  and  weighed.  Then  each  fragment  was  placed  in  a  device  that  measures  the  presented  area 
from  16  different  orientation  angles.+  Knowing  the  mass  of  the  fragment  and  its  presented  area  allows  us  to  derive  a 
value  of  the  shape  factor  by  making  use  of  equation  (1). 

6.2  Analysis  of  Shape  Factor  Data 

Figure  15  shows  the  shape  factor  distribution  for  one  of  the  spall  fragments,  as  measured  from  16  viewing 

angles. 


The  firings  and  fragment  collection  were  performed  by  Robert  F.  Kinsler  at  the  U.S.  Army  Research-Laboratory  Experimental  facility.  The 
measurements  of  the  fragment  presented  area  from  1 6  different  viewpoints  were  performed  at  the  U.S.  Army  Aberdeen  Test  Center. 
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Figure  15.  Shape  Factor  of  a  Spall  Fragment  from  16  Viewing  Directions. 

Although  the  frequency  distribution  shows  scatter,  the  shape  factor  is  relatively  uniform  over  orientation  angle. 
Since  16  viewpoints  are  not  very  many,  we  should  expect  a  certain  amount  of  variation.  For  example,  Figure  16 
shows  the  shape  factor  simulation  results  from  viewing  a  randomly  oriented  cylinder  with  an  LID  =  0. 1. 

3 

2 
N 

1 

0 

12  3  4 

r 

Figure  16.  Simulation  of  Shape  Factor  Distribution  from  16  Viewing  Directions. 


In  this  particular  shot,  150  spall  fragments  were  recovered.  They  ranged  in  mass  from  0.02  grains  to  several  hundred 
grains  (where  the  larger  masses  may  very  well  be  broken  pieces  of  the  penetrator).  The  fragment  masses  appear  to 
be  lognormally  distributed.1  Now  we  might  expect  there  to  be  a  correlation  between  the  mass  of  a  fragment  and  its 
shape  factor  distribution.  For  example,  one  might  expect  a  smaller  fragment  to  be  more  compact,  like  a  sphere,  and 
a  larger  fragment  to  be  more  irregular  in  shape.  The  first  question  to  address,  then,  is  whether  there  is  a  systematic 
trend  between  the  mass  of  the  fragment  and  its  shape.  To  test  this,  we  compute  the  linear  correlation  coefficient. 
The  correlation  coefficient  between  two  variables  x  and  y  is  defined  by 


cov  (x,  y) 

OxOy 


(42) 


t  The  behind-armor  debris  mass  distribution  is  commonly  fitted  with  a  Weibull  distribution,  where  ft  goes  by  the  name  of  Rosin-Rammler 
distribution.  We  have  found,  however,  that  a  lognormal  gives  a  better  fit  to  this  data. 
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where 


1  « 

cov  (x,  y )  = - r  Xte  "  *)(?/  “  y) 

n  1  I  —  i 


(43) 


is  the  sample  covariance  of  the  n  (x,  y)  data  pairs  and 

(44> 

n  1  i  =  i 

is  the  sample  variance.  The  value  for  r  always  lies  in  the  range  [-1, 1],  with  a  value  of  0  indicating  the  variables  are 
uncorrelated.  For  a  nonzero  value  close  to  zero,  we  use  the  fact  that  the  test  statistic 


Hn-2 
Vl  -r2 


(45) 


is  distributed  according  to  Student’s  distribution  and  use  it  to  test  the  null  hypothesis  that  the  actual  correlation  coef- 
ficient  is  zero. 


H0:p  =  0, 


(46) 


where  p  is  the  actual  correlation  coefficient.  The  computed  correlation  coefficient  between  the  mass  of  the  fragment 
and  the  average  shape  factor  is  found  to  be  r  =  -0. 130.  Our  test  statistic  is  t  =  - 1  •  60.  From  a  table  of  Student’s 
distribution,  we  find  that  this  is  smaller  (in  absolute  value)  than  1.96,  the  critical  value  for  a  two-sided  test  at  the  0.05 
'level  of  significance.  This  means  we  cannot  reject  the  hypothesis  that  these  two  variables  are  uncorrelated.  Another 
indicator  of  compactness,  or  lack  thereof,  is  the  difference  between  the  maximum  and  minimum  shape  factors.  The 
correlation  coefficient  between  this  difference  and  mass  is  found  to  be  r  =  -0. 076.  The  test  statistic  is  t  -  -0.93. 
Once  again,  since  this  is  smaller  than  1.96  in  absolute  value,  we  cannot  reject  the  hypothesis  that  these  two  variables 
are  uncorrelated.  These  tests  show  that  there  does  not  appear  to  be  any  statistical  correlation  between  shape  factor 
and  mass.  Therefore,  we  do  not  introduce  any  systematic  bias  by  pooling  all  the  shape  factors  regardless  of  mass. 
This  gives  us  a  total  of  2,400  shape  factors  (i.e.,  16  views  x  150  masses)  from  which  to  determine  the  shape  factor 
distribution.  This  is  shown  plotted  as  a  histogram  in  Figure  17. 


Figure  17.  Empirical  Shape  Factor  Probability  Distribution. 


The  empirical  distribution  appears  to  be  lognormal — and  resembles  very  much  the  distribution  obtained  by  varying 
the  LID  ratio  of  disk-like  cylinders  (cf.  Figure  13). 
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6.3  Monte  Carlo  Simulation  of  Shape  Factor 

Thus,  it  appears  that  randomly  oriented,  disk-like  cylinders  have  the  potential  for  simulating  real  behind- 
armor  debris  fragments.  This  is  about  as  far  as  we  can  go  without  additional  experiments.  The  experiments  cited  in 
this  report  do  not  attempt  to  classify  fragments  as  either  spall,  originating  from  the  RHA  target,  or  broken  pieces  of 
the  penetrator.  The  penetrator  in  this  case  was  made  of  tungsten  with  a  material  density  of  17.6  gm/cm3.  If 
denotes  the  fragment  shape  factor  computed  assuming  that  the  material  was  steel,  and  the  fragment  was  actually 
composed  of  tungsten,  then  the  correct  shape  factor  is 

nv  =  ySt(17.6/7.8)2/3  =  1.72ySt.  (47) 

So  this  can  have  a  significant  effect  upon  the  shape  factor  distribution.  Indeed,  at  face  value,  we  cannot  obtain  an 
LID  ratio  for  all  of  the  empirical  data,  due  to  the  fact  that  the  smallest  average  shape  factor  for  a  cylinder  is  1.3895. 
The  smallest  average  shape  factor  from  the  empirical  data  is  0.8424.  Notice,  however,  that  if  this  is  actually  a  tung¬ 
sten  fragment,  then  the  average  shape  factor  is  1 . 72  x  0. 8424  =  1 . 4489,  which  does  lie  within  the  range  of  a  cylin¬ 
drical  shape. 

A  procedure  for  using  disk-like  cylinders  to  simulate  behind-armor  debris  fragments  can  be  formulated  as 

follows: 

•  Collect  behind-armor  debris  fragments  and  identify  the  material  density. 

•  Use  equation  (15)  to  compute  the  LID  ratio  from  the  average  shape  factor. 

•  Identify  a  probability  density  function  that  captures  the  LID  distribution.  Denote  this  distribution  by 
F{LID)\  the  simulation  leading  to  Figure  13  assumed  that  this  was  the  Uniform  distribution. 

Once  the  F(LID)  distribution  is  known,  an  algorithm  for  a  Monte  Carlo  simulation  is  as  follows: 

(1)  Generate  LID  ~  F(LID). 

(2)  Generate  (0,  <p)  ~  Surface  of  the  Unit  Sphere. 

(3)  Return  y  =  a  sin  0  4-  b  cos  0,  where  a  and  b  are  given  by  equation  (7). 

Once  the  mass  m,  density  p,  and  LID  ratio  of  the  fragment  are  known,  then  the  actual  dimensions  of  the  fragment 
can  be  computed: 

DJ  m/p  T  and  L  =  (L\d.  (48) 

U*/4)(LAD)J  U>J 


7.  CONCLUSIONS 

A  dimensionless  shape  factor  was  defined  and  applied  to  a  randomly  oriented  right  circular  cylinder.  This 
analysis  established  the  following  results. 

•  The  LID  ratio  of  a  disk-shaped  cylinder  can  be  obtained  from  the  average  shape  factor  (see  p.  6). 

•  An  explicit  probability  density  function  can  be  obtained  for  a  randomly  oriented  cylinder. 

•  Two  algorithms  can  be  used  for  simulating  randomly  oriented  cylinders: 

□  A  fast  method,  derived  from  the  inverse  transformation  technique. 

□  An  explicit  simulation  of  the  oriented  cylinder. 

•  Shape  factors  for  disk-shaped  cylinders  show  potential  for  simulating  behind-armor  debris  fragments. 

It  should  also  be  noted  that  the  average  shape  factor  should  not  be  used  in  place  of  the  shape  factor  distribu¬ 
tion.  For  example,  consider  a  1-gm  steel  cylinder  with  an  LID  =  0. 1  striking  a  1/32-inch  mild  steel  plate  at  a  speed 
of  360  m/s.  The  THOR  penetration  equations  indicate  that  the  projectile  will  perforate  the  plate  only  if  y  <  2. 3. 
Since  the  average  shape  factor  for  this  LID  ratio  is  2. 6,  using  y  to  represent  the  projectile,  will  lead  to  the  conclu¬ 
sion  that  there  is  no  perforation.  Using  the  shape  factor  distribution,  however,  shows  that  there  is  indeed  perforation 
in  42%  of  the  impact  encounters  due  to  the  effect  of  random  orientation  (cf.  Figure  12). 

The  results  derived  in  this  report  should  be  used  for  designing  fragment-simulating  projectiles.  And  since 
the  FATEPEN  code  explicitly  models  penetrators  as  right  circular  cylinders,  fully  accounting  for  both  orientation 
and  LID  ratio,  these  results  should  also  be  useful  when  using  that  model  for  simulating  fragment  penetration. 
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APPENDIX  A: 

MEAN  PRESENTED  AREA  OF  A  CONVEX  SOLID 


Let  n  be  the  outward  normal  on  the  surface  of  the  solid  and  let  n'  be  a  unit  vector  along  an  arbitrary  direc¬ 
tion.  The  presented  area  of  the  convex  solid  in  the  direction  n'  is  given  by  the  integral 


Ap=  J  n  •  iT  dS, 

n-n'  >  0 


(A-l) 


where  dS  is  an  element  of  surface  area.  The  mean  presented  area,  averaged  over  all  directions  iT,  is  given  by 

Ap  =  ~j  dQ'  J  n  n'dS,  (A-2) 

n  n'  >  0 

where  the  first  integral  is  over  all  solid  angles  of  a  unit  sphere.  Interchanging  the  order  of  integration,  we  have 

Ap  =  jdS-^  J  n-n'dQ'.  (A-3) 

n  n'  >  0 

In  the  last  integral,  n  remains  fixed  while  n'  varies  over  all  solid  angles  such  that  n  •  n'  >  0.  By  changing  coordinates 
so  that  n  is  along  the  polar  axis,  n  •  it'  =  cos  0 ,  and  this  integral  is  easily  evaluated: 

2  ix  n!2 

J  n-n'  dQf  =  J  d<p  J  cos  0  sin  OdO  =  n .  (A-4) 


n-n'  >0  0  0 

Substituting  this  result  into  equation  (A-3)  gives 


(A-5) 


The  integral  is  the  total  surface  area  of  the  solid,  so  we  have  the  general  result: 


,t 


The  mean  presented  area  of  a  convex  solid  is  one-fourth  of  the  total  surface  area. 


It  is  a  well-known  fact  that  a  sphere  is  the  shape  that  minimizes  surface  area.  With  the  aid  of  the  above 
result,  we  can  conclude  that  the  minimum  mean  presented  area  of  any  convex  solid  is 


1  2 
Ap  >.  -  ( surface  area  of  a  sphere )  =  nr\ 


(A-6) 


where  r  is  the  sphere  radius.  Substituting  this  into  the  equation  Ap  =  /(m/p)273,  we  find  that  the  average  shape  fac¬ 
tor  for  any  convex  solid  obeys  the  inequality 


7  >  (3/2)2/3(tt/4) 1/3  «  1.20899. 


Notice,  in  particular,  that  a  right  circular  cylinder  obeys  this  inequality  since  it  has  a  minimum  average  shape  factor 
of  (3/2)(;r/4)1/3  »  1.38395. 


This  theorem  has  been  attributed  to  Cauchy.  It  was  first  brought  to  the  author’s  attention  by  Dr.  Robert  Shnidman  who  also  outlined  this 
particular  derivation. 
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APPENDIX  B: 

SHAPE  FACTOR  OF  A  SPIN-STABLIZED  FRAGMENT  SIMULATOR 

The  dimensions  of  the  chisel-nose  fragment  simulator1  are  shown  in  Figure  B-l. 


D/3 


D/15 


Top  View 

Figure  B-l.  Dimensions  of  Fragment  Simulator. 


Shape  Factor 

The  fragment  is  spin-stablized  so  that  the  chisel  nose  is  forward.  The  dimensionless  shape  factor,  y,  is 
defined  by 

Ap  =  y(m/p)2'3,  (B-l) 

where 

Ap  is  fragment  presented  area  (cm2), 
m  is  fragment  mass  (gm),  and 

p  is  fragment  density  (gm/cm3). 

Volume 

The  fragment  volume,  mlp ,  can  be  calculated  as  follows: 

Total  Volume  =  Volume  of  Cylinder  with  Diameter  D  and  Height  2 D/5  -  Volume  of  Chiseled  Edges.  (B-2) 
The  volume  of  the  chiseled  edges  can  be  calculated  using  the  coordinate  system  shown  in  Figure  B-2. 


1  Fisher,  Todd  J.  and  Robert  Shnidman.  “An  Empirical  Model  for  Predicting  Electrical  Cable  Failures  in  Modem  Combat  Vehicles,”  ARL- 
MR-28,  U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  January  1993. 
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Figure  B-2.  Coordinate  System  for  Calculating  Volume  of  Chiseled  Edges, 


Carrying  out  the  integrations,  we  find 

m  \k  19V2  sin~1(2V2/3)l  3 

p  ~  10  162  +  12 


Substituting  this  into  equation  (B-l),  with  Ap  =  rtD2!A,  gives 


7  =  T 


it  ["  k  19V2  sin  1  (2V2/3) 


4  10  162 


«  1.97448. 


(B-3) 


(B-4) 


(B-5) 


Shape  Factor  of  a  Spin-Stablized  Right  Circular  Cylinder 

The  shape  factor  of  a  (flat  nose)  right  circular  cylinder  that  impacts  face-on  is 


Thus,  the  proper  LID  ratio  to  effect  a  given  shape  factor  y,  is 

LID  =  V^47"3/2. 


(B-6) 

(B-7) 
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City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the  Old 
or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 
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